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ABSTRACT 

We present a statistical comparison of the mass ratio distribution of companions, as observed in 
different multiplicity surveys, to the most recent estimate of the single object mass function (Bochanski 
et al. 2010). The main goal of our analysis is to test whether or not the observed companion mass 
ratio distribution (CMRD) as a function of primary star mass and star formation environment is 
consistent with having been drawn from the field star IMF. We consider samples of companions for 
M dwarfs, solar type and intermediate mass stars, both in the field as well as clusters or associations, 
and compare them with populations of binaries generated by random pairing from the assumed IMF 
for a fixed primary mass. With regard to the field we can reject the hypothesis that the CMRD was 
drawn from the IMF for different primary mass ranges: the observed CMRDs show a larger number 
of equal-mass systems than predicted by the IMF. This is in agreement with fragmentation theories 
of binary formation. For the open clusters a Persei and the Pleiades we also reject the IMF random- 
pairing hypothesis. Concerning young star-forming regions, currently we can rule out a connection 
between the CMRD and the field IMF in Taurus but not in Chamaeleon I. Larger and different samples 
are needed to better constrain the result as a function of the environment. We also consider other 
companion mass functions (CMF) and we compare them with observations. Moreover the CMRD 
both in the field and clusters or associations appears to be independent of separation in the range 
covered by the observations. Combining therefore the CMRDs of M (1-2400 AU) and G (28-1590 
AU) primaries in the field and intermediate mass primary binaries in Sco 0B2 (29-1612 AU) for mass 
ratios, q = M2/M1, from 0.2 to 1, we find that the best chi-square fit follows a power law dN/dq oc , 
with /? = —0.50 ±0.29, consistent with previous results. Finally we note that the KS test gives a ~1% 
probability of the observed CMRD in the Pleiades and Taurus being consistent with that observed 
for solar type primaries in the field over comparable primary mass range. This highlights the value of 
using CMRDs to understand which star formation events contribute most to the field. 

Subject headings: binaries: general - Stars: formation - Stars: pre-main sequence - open clusters and 
associations: general - Stars: luminosity function, mass function 



1. INTRODUCTION 



The study of binary stellar systems and their proper- 
ties is one of the most important topics in star formation. 
Since mai iy stars form in mult iple syst ems, both in the 
field (e.g iDuauennov fc Mavoif Il99ll: iFischer k, Marcv 
1992) and in star-forming regions (e.g iPatience et al 



20021 ) ■ a correct understanding of the binary fraction 



and companion mass distribution for different primary 
masses represents a fundamental test for star forma- 
tion theories. Binary populations, in fact, may carry an 
even wider amount of informa tion concerning star forma- 
tion p rocesses than the IMF ([Goodwin fc Kouwenhovenl 
l2009f ). By conventional definition, in a binary system of 
stars with masses Mi and M2, with Mi > M2, Mi and 
M2 indicate the primary and secondary mass, respec- 
tively. Consequently one can define the ratio of the sec- 
ondary over the primary mass as q — M2 /Mi . Analogous 
to the initial mass function (IMF) for single objects, one 
can define the companion mass ratio distribution (here- 
after CMRD) as the distribution of q for a chosen primary 
mass. 

Different binary formation models predict different 
mass ratio distributions and dependencies of the CMRD 



on the primary mass. Traditionally these classes of 
models have been divided into capture and fragmen- 
tation scenarios. Capture refers to the tidal capture 
of two unbound objects on a timescale that is long 
compared to th e collapse time of each component (e.g. 
iMcDonald fc C larke 1993). For each primary star the 
mass of the secondary is chosen randomly from the sin- 
gle star mass function and the secondary-mass distri- 
bution would reflect the IMF. While tidal capture ap- 
pears to be too inefficient in reproducing high binary 
fractions, it has been noticed that, particularly in small 
groups of stars , star-disk encounters may form binaries 
(jMcDonald fc C larke 1995). In any case even this disk 
assisted capture, whereby a star passing through the disk 
of another which dissipates enough kinetic energy to form 
a bound system, is unlikely to be the most re levant bi- 
nary formation mechanism (iBoffin et al.lll998l ). 

Fragmentation scenarios are the preferred mechanism 
for the formation of multiple systems. The so-called 
fragmentation models are u s ually cl assified a s prompt 
fragm entation (e.g. iBossI Il986t iBonnell fc Bastie^ 
19921) and disk fragme i itatio n (e.e. iBonnelL .199 
Stamatellos fc WhitworthI [2OO90 . In the prompt frag- 
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mentation scenario both primary and secondary stars 
form by fragmentation of the same collapsing molecular 
cloud core. Disk fragmentation takes place in a newly 
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formed star-disk system in which the disk subsequently 
fragments due to density perturbations. The latter mech- 
anism is a process by which low-mass stars and BDs (e.g. 
companions with q < 0.25) may form. In both cases frag- 
mentation is only the first step in binary formation and 
processes such as disk accretion and dynamical interac- 
tions all contri bute to dete rmine the final properties of 
binary systems ()Batell200l) . In general, continued accre- 
tion onto both objects from a common reservoir tends in 
the long term to equalize the masses, moving the q dis- 
tribution towards unity, and this effect seems to be more 
significant f or high mass primaries and in closer binaries 
(|Batdl2000n . 

Recently, a variation on capture has been pro- 
posed as mechani sr n for forming wide b inarie s 
(jKouwenhoven et all [20Tot iMoeckel fc Rati [20T0I ). 
In this scenario wide binaries (10'' - 10^ AU) would 
form during the dissolution phase of star clusters, 
especially during the quick expansion of clusters after 
gas expulsion. This mechanism could perhaps explain 
the substantial population of wide binaries observed in 
the field and, to the first order approximation, the mass 
ratio distribution could be similar t o that expected from 
rando m pairing of individual stars (jKouwenhoven et al.l 

However, at the moment no model is able to re- 
produce all of the observed binary properties, in par- 
ticular, the predicted distributions of separations and 
mass ratios tend not to m atch the observations very 
well ([Goodwin et al.l 120071 ). Even if no observation 
will definitely confirm one theory, observations of CM- 
RDs in disparate environments can at least put con- 
straints on theoretical binary formation models. Begin- 
ning with the pioneering multipl icity survey of G s tars i n 
the solar neighborhood by Duquennov fc Mavoii ()199H ). 
there have been several studies o f binar y properties in 
the fi eld (e.g. iFischer fc M arcV ('1992'), 'Reid fc Gizi^ 
(I1997D for M dwarfs. iMetche v & Hillenbrand (200^) and 
iRaghavan et al.l (|2010l l for solar- tvpc stars) and i n clus- 
ters or associations (e.g. [ L eincrt ct al. 1993; Ghc z et al.l 
ll993tlPatience et al.ll2002HKraus et al.ll201lD . Curiously 
in young clusters the binary fractions overall are higher 
but little is known about the CMRD. 

Furthermore observations of different star-forming re- 
gions revealed that the mass dis tribution of cores usually 
has a form simila r to the IMF (Motte et al.lll998l 120011: 
lAlves et aLll2007D . leading to the suggestion that IMF 
and the core mass function are directly related. At the 
same time the majority of observations of single objects 
from the field, local young clusters, old gl obular clusters 
and a ssociations suggest a universal IMF (jBastian et al.l 
I2010f) . What is the role played by binaries and multi- 
ple systems? Is there a connection be t ween the CMRD 
and the IMF? IMetchev fc Hillenbrand l (120091) present an 
analytical form of the companion mass function (CMF) 
and reject the hypothesis that the CMRD of solar-mass 
stars in the field is consistent with having been drawn 
from random pairing of the single star IMF. However, 
the most rece nt estimate of the log- normal IMF for iso- 
lated objects (jBochanski et al.|[2010[) is peaked at higher 
masses than those they considered, and might be in closer 
agreement with the CMRD they derived. It is also inter- 
esting to study CMRDs as a function of primary mass as 



well as in different environments and look for variations 
in the CMRD among regions where the IMF varies. 

Finally, another crucial question that the CMRD could 
address is the origin of the field. Matching properties of 
the field to star-forming regions could give insights into 
what sort of regions contribute most to the field. Bi- 
naries are subject t o dynamical evolution and disruption 
(jParker et al.lf2009l ). changing the overall binary fraction. 
The extent of these dynamics depends on the environ- 
ment in which they were born. Recent results from N- 
body simulations also predict that some binary proper- 
ties (e.g. the CMRD for low-mass bi naries) might be in- 
depen dent of dynamical processing (jParker fc Goodwinl 
l2011f ). making them excellent tracers of origins. 

A rigorous approach to answering these questions re- 
quires a careful account of the completeness of obser- 
vational data and potential biases as a function of sep- 
aration, besides a proper choice of the IMF. So far, a 
complete analysis of the CMRD over a broad range of 
primary masses and as a function of separation and en- 
vironment has not been done. In this paper we address 
the problem of the connection between IMF and CMRD 
by considering samples of binaries with primaries of dif- 
ferent masses and using the most recent evaluation of 
the IMF. We first describe the datasets we have used in 
our analysis (Section [5]). Then in Section [3] we discuss 
the methodology we adopt while the results obtained are 
shown in Section H) Finally, Section [5] and Section [6] are 
left to the discussion of the results and to our conclusions. 

2. DATASETS 
2.1. Samples from the field and Sco 0B2 

As mentioned before, the multiplicity of stars in the 
field has been investigated in the past years by different 
groups. Among these surveys we have selected for our 
analysis three studies, each of which surveyed compan- 
ions for a restricted range of primary masses (M dwarfs, 
solar- type stars and intermediate mass stars). In this 
section we give a brief overview of the datasets we have 
chosen while in section 12.2! we describe the samples of 
binary systems from young clusters/associations that we 
have also considered in our analysis. A full summary 
of the main properties of all these samples is given in 
Table [1] We already said that to obtain reliable results 
it is important to account for completeness and possi- 
ble biases. However also the estimate of completeness 
and observational biases is not free from uncertainties. 
Therefore, instead of considering completeness-corrected 
samples, we decided to limit our investigation to the mass 
and separation range where the completeness of the sam- 
ples is flat (so the sample is representative in mass, if not 
complete) and exceeds a certain level (> 65%). 

The sample of M dwarfs we have considered is t he set 
of binary systems collected by IFischer fc Marcvl ()1992f ) 
(hereafter FM92) from several high-quality surveys of M 
dwarfs with distances within 20 pc from the Sun (ages 
> Gyr) . Each one of these surveys covers a different an- 
gular separation range, but the complete sample extends 
from roughly 1 to 2400 AU in separation and down to 
q = 0.2 in mass ratio. Generally, M dwarfs with masses 
< 0.2 Mq show mass r atios biased towards un ity due to 
sensitivity limitations (I Fischer fc Marcvlll992D . For this 
reason we have considered only binary systems with pri- 
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TABLE 1 

Sample Properties 



Sample Rcf ^ Primary Type No. Multiple Systems Separation Range (AU) qu. 



Field 1 M 27 1-2400 0.2 

Field 2 F/G 30 28-1590 0.1 

ScoOB2 3 A/late-B 60 29-1612 0.05 

Pleiades 4 F/G 22 11-910 0.2 

a Persei 5 F/G 18 26-581 0.25 

Chamaeleon I 6 G/K^ 13 20-800 0.1 

Taurus 7 G/K" 40 5-5000 0.1 



References: (1 ) IFischer fc Marcvl (TM^) . (2) IMetchev fc Hillenbrandl f200g), (3) 
IKouwenhovcn ct al. (20051 ), (4) IBouvier et al.l OWm . (5) IPatience et al.l I I2002I ). (6) 
[Talrcnicre et al. (2008), (7) IKraus eFall i20TH ). 
^ The mass range is 0.55 and 2.2 Mq, comparable to MH09 
The mass range is 0.7 and 2.7 Mq, comparable to MH09 



maries having masses between 0.2 to 0.55 Mq, where the 
sample is 85% complete. The sample consists then of 27 
systems. 

Regarding solar-mass stars we s e lected the work 
presented in IMetchev fc Hillenbrandl ()2009f ) (hereafter 
MH09). MH09 report results from an adaptive op- 
tics survey of stellar and substellar companions to solar 
analogs (range in primary mass between 0.7 to 1.3 Mq) 
within 10-190 pc and in the 3 Myr - 3 Gyr age range. The 
orbital separation interval covered is 28-1590 AU. The 
choice of this survey, with r espect to previous works (e.g. 
IDuauennov fc Mavoiill99l[) , is due to the higher sensitiv- 
ity to low-mass companions, meaning small mass ratios 
{q < 0.1). In order to have a 65% complete sample we 
considered the set of 30 binary systems with q > O.l and 
companions between 28 and 1590 AU from the primary 
which was defined as minimally-biased sample (AD^q) in 
their paper. 

Finally, we chose a sample of companions to A-type 
and late B-type primaries. Due to the shorter lifetime 
of more massive stars and the difficulty to find a statis- 
tically large and complete survey of intermediate mass 
primaries in the field, we selected a dataset in the young 
(5-20 Myr) and nearby (~140 pc) Scorpius 0B2 associ- 
ation (ScoOB2). This binary population was observed 
in the near-infrared adaptive optics multiplicity survey 
described by Kouwenliovcn ct al. (2005,) (hereafter K05) 
and the properties of the 60 stell ar systems we have 
used in our analysis are taken from iKouwenhoven et ahl 
([2007). This survey is sensitive to very low mass ratios 
(down to g ~ 0.05) over the orbital separations range 
29-1612 AU between primaries and companions. Despite 
the fact that Sco 0B2 is not a sample from the field and 
it is young, it is still the best sample (> 90% complete) of 
intermediate mass primary binaries we can study. There- 
fore we will include it in the analysis of the other datasets 
from the field. 

We emphasize that the three datasets span similar sep- 
aration ranges. Even though the study of the detailed 
dependence of the CMRD on angular separation goes be- 
yond the purpose of the present work, in section [4.11 we 
will show that our results should not be affected by any 
possible change in the shape of the CMRD with orbital 
separation. 



We considered also four datasets of companions to 
solar-type stars from nearby clusters or associations that 
over similar separation ranges have a reasonable number 
of binary systems. We have selected observations of bina- 
ries in two open clusters, Pleiades and a Persei (a Per), 
and in two T associations, Chamaeleon I and Taurus. 

The Pleiades is one of the best studied open clusters, 
due to its proximity and richness (roughly 1000 stars 
at a distance of ~ 120 pc). With an ag e of 125-150 
Myr (iStauffer et al.lll998tlBiIrke et al.ll2004l) it is just old 
enough t o be dynamically evo lved. The sample from the 
Pleiades (jBouvier et al.l[l997l ) consists of 22 binary sys- 
tems with G and K primaries observed in the near-IR 
using adaptive optics. The separation range covered by 
this survey is 11-910 AU and the mass ratio distribu- 
tion is more than 70% complete down to 0.2 over this 
separation range. 

a Per is a younger c luster, with an age of ~90 
Myr IS tauffer et all I1999D . at a distance of ~ 190 pc 
(jRobic hon ct al.l ll999l ). We selected a sa niple of 18 solar 
type stars within the dataset presented in IPatience et all 
(2002). They were nearly complete in the separation 
range from 26 to 581 AU and were sensitive to mass ra- 
tios q > 0.25. 

C hamaeleon I, i nstead , is one of the nearest (^170 
pc, Bcrtout et al.l Il999t ) low-density young (^1 Myr, 
iLuhman 2004) star-forming regions. It consists of ^230 
stars and has a stella r density that is low compared to 
other young regions (jLuhmanl I2008D . We have consid- 
ered the results of a multiplicity survey presented by 
iLafreniere efall (|200l ). The primaries span the mass 
range from 0.1 to 3 Mq and the separation range ~20- 
800 AU. We have selected a subsample with only K and 
G primary binaries with masses between 0.55 and 2.2 
M© and mass ratios down to g 0.1 (^ 90% complete), 
comparable to MII09 (13 systems in total). 

Finally we selected a sample of solar-type primary bi- 
naries in Taurus f rom the almost complete sample by 
iKraus et al.l (|201lD . Taurus is another young (1 Myr) 
low-density star-forming region close to the Sun (d — 
140 pc) with niore th an 300 PMS stars and brown dwarfs 
(jKenyon et al.ll200"8l ). We considered 40 systems with pri- 
mary masses between 0.7 and 2.5 Mq, mass ratio q > 0.1 
and angular separation in the range 5 — 5000 AU. 

We will discuss the dependence of the CMRD on sep- 
aration in section 14.11 



2.2. Clusters or associations 
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3. METHODOLOGY FOR OUR ANALYSIS OF THESE 
SURVEYS 

3.L Monte Carlo simulations 

Our goal is to explore whether the observed mass ratios 
in the field and young clusters or associations as function 
of primary mass could be the outcome of random pairing 
of stars from the stellar IMF. To this end we have created 
a Monte Carlo tool able to generate artificial companion 
mass ratio distributions as expected by random sampling 
of secondaries from a chosen function, for fixed primary 
mass. Through these simulations one can reproduce a 
population of N binaries by fixing the mass of the primary 
and the analytic form to be tested, and then compare this 
simulated CMRD with the observations. 

3.2. Initial Mass Function 

The Initial Mass Function we have c onsidered in our 
analys is is the single objects IMF from IBochanski et al.l 
(|2010f l (hereafter Bo2010). Below 1 Mq it is a log- 
normal function of the mass, defined as ^(logm) — 
dn/d\ogm. Except for a normalization constant, it can 
be parametrized (in (log Mo)^^pc^'^) as: 

^ / (l0gTO-l0gTOc)^\ ... 

C(log m) oc exp j ^ '> (1) 

where nic = 0.18 and a = 0.34. 

For m>l M© we assumed the classical "Salpeter 
slope" : 

^(logm) oc m"^-^^. (2) 

The study by IBochanski et al.l (120101) is ba, s ed on 
the observational work presented in ICovev et al.l (|2008[ ) 
which represents the largest field investigation of the lu- 
minosity function to date constructed from a catalog of 
matcheii Sloan Digital Sky Survey (SDSS) and Two Mi- 
cron All Sky Survey (2MASS) sources. Note that in 
Bo2010 the log-normal peak of the mass di stribution 
is shif ted toward higher masses compared to iChabrieil 
(pool (mc = 0.08 and a = 0.69). 

In each run of our Monte Carlo simulations, the as- 
sumed IMF is normalized to the primary mass that we 
choose, to the appropriate range of q for the dataset with 
which we compare the results, and to the number of bi- 
naries N that we want to reproduce. We typically run 
each simulation 10^ times. 

3.3. KS test 

To evaluate the probability that the observed CMRDs 
and the simulated ones come from the same parent distri- 
bution we apply the Kolmogorov-Smirnov test (KS test) 
to the cumulative distributions of q values. The KS test is 
a statistical test which returns the probability that two 
distributions were drawn from the same parent sample 
by examining the maximum difference in the cumulative 
distribution functions. 

For our purposes, we have tested how well this statis- 
tical tool can distinguish differences in the shape of the 
distributions and evaluated the extent to which results 
depend on sample size. 

First, we tested the reliability of the KS test to dis- 
tinguish two populations of stars distributed in mass ac- 
cording to power law distributions with different slopes 
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Fig. 1. — Capability of the KS test to distinguish flat and log- 
normal distribution. The solid black line describes how varies 
the probability given by the KS that a log-normal differs from 
a flat distribution, as function of sample size. We computed 
this probability for samples of 10, 20, 30, 40, 50, 60, 70 and 
80 objects. 

as a function of sample size. With larger populations, 
the KS test is able to detect smaller differences in slope. 
From the comparison of distributions of only 10 objects 
the KS test gives a probability of ~ 10~^ for a differ- 
ence of 5 in the slopes while when the number of objects 
increases to 30 the same probability is already obtained 
with a difference of 2.5. 

Second, we checked the capability of the KS test to 
distinguish a log-normal and a flat distribution, again as 
function of sample size. This test is of great importance 
because, on one hand, we want to test the hypothesis of a 
log-normal CMRD, on the other, the linearly flat distri- 
bution of g is a c ommonly made assumption in numerical 
simulations (e.g. iKouwenhoven et aLll2009t iParker et al.l 
l2009h . In Figure ID we show our results. In this compar- 
ison, with a sample of ~50 objects the KS test returns 
a 1% chance of having been drawn from the same par- 
ent. A KS probability of 1% is the threshold we adopt 
equal to or below which we reject the hypothesis of two 
distributions being consistent. 

4. RESULTS 

Here, we report our flndings regarding the comparison 
of the observed CMRDs with the simulations. We begin 
describing in section H?T1 the results for the samples of M 
dwarfs, G stars in the held and intermediate mass stars 
in ScoOB2. In section 14.21 we summarize the outcome 
of our tests for the Pleiades, a Per, Chamaeleon I and 
Taurus. We compare the observed CMRDs with other 
commonly assumed companion mass functions in section 
14.31 Finally (section 14. 4p we give our best-flt estimate 
of the combined distribution of M dwarfs and G stars in 
the field and intermediate mass stars in ScoOB2. 

4.1. Results from the field and Sco 0B2 

The top panel of Figl2] shows the CMRD for the 
sample of 27 M dwarf primary binary systems from 
iFischer fc Marcvl (|l992). The hatched histo gram rep- 
resents the observed distribution of q, while the dashed 
line is the CMRD generated by random pairing through 
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Monte Carlo simulations from Bo2010 for the same range 
of mass ratios {q >0.2). The KS test gives a probability 
of ^ 1% that the observations are consistent with the 
IMF in the separation range 1-2400 AU. From the fig- 
ure it appears that there is an overabundance of equal 
mass binaries in the observed sample compared to the 
predictions of random pairing from the IMF. 

In Figl2] (central and bottom panels) we also present 
the results obtained for the other two datasets. In both 
cases we compare the observations (hatched histograms) 
with the simulated CMRDs (dashed line histograms) over 
the same range of q. On the basis of the KS test for the 
samples of solar type from the field and A type stars from 
Sco 0B2 we get probability of 10~^ and 10"^'^ respec- 
tively that the CMRD is consistent with the field IMF 
in the separation range ~30-1600 AU. The results are 
summarized in Table [2l We find again that the observed 
mass ratios are more strongly peaked towards unity than 
in the simulations from the IMF. 

This overall result suggests that we can reject the hy- 
pothesis that the CMRD is consistent with having been 
drawn from the IMF over the separation range ^^30-1600 
AU and this statement seems to hold true independent 
of the primary mass and angular separation. In fact we 
checked with the KS test whether for each sample we 
see variations in the CMRD as a function of the angular 
separation. Practically, for each one of the datasets we 
have considered different values of the angular separa- 
tion within the range covered by the observations, and 
we evaluated for each of these separations the KS test 
probability of the CMRD inside this value being consis- 
tent with the CMRD outside. For any given separation 
we find probabilities less than 0.1%. Therefore we did not 
see any evidence for dependence on orbital separation in 
any of the samples under study. 

We also note that the CMRDs generated for different 
primary masses through Monte Carlo simulations and 
shown in Fig. [3] are significantly different. This is ex- 
pected since the random pairing from the IMF predict a 
strong dependence of the CMRD on the primary mass. 
For example, a CMRD drawn from the IMF for primary 
mass near 1 Mq should exhibit a peak near q—0.18, 
whereas for primary masses near 0.2 Mq should decrease 
monotonically, as shown in Fig. (3] 

As we do not see variation of the CMRD with angu- 
lar separation, we can compare with the KS test also 
the observed CMRDs for these three datasets over the 
common range of q. The distributions of mass ratios for 
M dwarfs and solar type stars are consistent at the 6% 
level. On the other hand the probability of the sample of 
intermediate mass stars from ScoOB2 being consistent 
with G/K stars and M dwarfs in the field is 36% and 
53% respectively. These results suggest that we cannot 
reject the hypothesis that they are all drawn from the 
same parent distribution. 

4.2. Results from clusters and associations 

We also compared the datasets of solar type stars from 
the Pleiades, a Per, Chamaeleon I and Taurus with the 
results of the MC simulations of the IMF over the range 
of q spanned by each sample. In Figure [3] we show the 
CMRD for the Pleiades. The KS test probabilities that 
the observed CMRDs of all samples are consistent with 
having been drawn from the IMF are given in Table [H 
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Fig. 2. — Companion mass ratio distributions and the IMF 
in the field. From top to bottom the comparison between the 
observed CMRD and the log-normal field IMF is shown for M 
dwarfs, G stars in the field and for the sample of intermediate 
mass stars in Sco 0B2 respectively. The hatched histogram 
represents the observed CMRD for the respective dataset of 
binary systems (see Section [5]). Superimposed with a dashed 
line is the CMRD generated for the same number of objects 
through random pairing from Bo2010. The KS probabilities 
are summarized in Table [2l 
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TABLE 2 

KS TEST PROBABILITIES 



Sample - Type Ref No. Systems Bo2010 (%) Flat CMF (%) dN/dM2 oc Mj""'" (%) 



Field - M 1 27 1 58 24 

Field - F/G 2 30 lO-^ 2 58 

ScoOB2 - A/late-B 3 60 IQ-^^ 0.4 30 

Pleiades - F/G 4 24 lO""* 34 17 

o Persei - F/G 5 18 0.1 27 89 

Chamaeleon I - K/G 6 13 17 30 76 

Taurus- K/G 7 40 10"" 45 2 



Referenc es: (11 IFischer fc Marcvl 1 11993) ■ (2')IMetchev fc Hillenbraiidl 1120091). fS) IKouwenhoven et al 
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Fig. 3. — Companion mass ratio distribution for solar-type 
stars in the Pleiades. The image shows the comparison be- 
tween the observed CMRD in the Pleiades and the predictions 
from the IMF. We adopt the same legend as in figured) The 
probability from the KS test that observations are consistent 
with the IMF is less than 1%. 

With regard to the Pleiades, a Per and Taurus (KS test 
probabilities of 10-''%, 0.1% and 10""%, respectively) 
we can reject the hypothesis of random pairing, while the 
higher probability (17%) between the IMF and CMRD 
in Chamaeleon I does not allow us to rule out the null 
hypothesis in this case. However, the number of objects 
(13) in Chamaeleon I sample is quite low. As we have 
shown in Section [3. 31 the KS test can only distinguish ex- 
treme differences in distributions from such small sam- 
ples. 

We should also note that in Taurus the IMF is peaked 
toward higher mass es with respect to the field IMF 
(jLuhman et al.|[2009l) . The use of the proper mass dis- 
tribution would bring the CMRD in Taurus in closer 
agreement with random pairing from the IMF. Hence 
we should be cautious in interpreting this preliminary 
result. 

To summarize, at the moment in the Pleiades and a 
Per we can reject the possibility of the CMRD being 
drawn from the IMF for orbital separation between 20 
and 600 AU, whereas concerning much younger regions, 
we rule out the random pairing from the field IMF only 
in Taurus in the separation range 5-5000 AU. Data from 
larger and different samples are needed to better con- 
strain the result as a function of age and environment. 



4.3. Different Companion Mass Functions 

Using the same Monte Carlo method, we have tested 
also whether the observed CMRD as function of primary 
mass and environment is consistent with other analytic 
forms of the CMRD. First of all we have considered a 
linearly flat companion mass ratio distribution (see Sec- 
tion [331) second we have tested the companion mass 
distribution dN/dA'h oc M2"°" suggested by MH09. In 
Table [2] we report the KS probabilities for each dataset. 

Concerning the flat distribution, only for the sample 
of A and late B-type primary binaries in Sco 0B2 we 
can reject the hypothesis that the two distributions are 
consistent. The comparison in this case is shown in the 
top panel of Figure H) Note that the ScoOB2 dataset is 
the largest sample, placing the strongest constraints on 
possible differences. The MC simulations of a flat CMF 
for the young regions match well the observations (see 
e.g. bottom panel of Figure |4]). Regarding the CMRD 
provided by MH09 we find a KS probability exceeding 
15% for all samples (see Table[2]) except for Taurus (2%). 

We should keep in mind that the KS test is not suited 
to evaluate which is the best fit distribution. If we take 
as an example the results for the FM92 sample, the dif- 
ference in the probability from 58% to 24% between the 
flat and M09 CMF in the context of the KS test does 
not have any significance. Furthermore, the sample size 
of our datasets in the majority of cases prevents us from 
discriminating between log-normal, flat or other distri- 
butions (see Section 13.31) . For this reason we have uti- 
lized a chi-square procedure to determine the best fit for 
the CMRD for a combined sample including all primary 
masses. 

4.4. Chi-square best fit 

Motivated by the fact that the CMRD appears to be 
independent of angular separation over the range we 
are considering and that the distributions are not dis- 
tinguishable, we combined together, over the common 
range of mass ratios (q=0.2-l), the samples of M dwarfs 
and G stars in the field and intermediate mass stars in 
ScoOB2 even though the separation ranges vary across 
the samples. We then used the composite q distribution 
to find the best fit. According to the chi-square test for 
Mi=0.25-6.5 Mq the total mass ratio distribution fol- 
lows a power-law dN/dq cx g^, with /3=-0.50±0.29 {x^— 
0.7 with 7 degrees of freedom; see Figure [S]). 

This result is also in agreement with the value 
of /3=-0.50 (45-900 AU) for B star primaries 
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Fig. 4. — Test of other CMFs Top: the figure shows the com- 
parison between the observed CMRD for intermediate mass 
stars in ScoOB2 with a flat CMRD (dotted line) and a com- 
panion mass function of the form dN/dM2 oc * (dashed 
hne). We found a probability of less then 1% that the obser- 
vations are consistent with the flat CMF while a 30% level of 
agreement with the CMF by ME[09. Bottom: Comparison for 
the observed CMRD for solar-type stars in the Pleiades with 
the two choices of CMF. The KS test probabilities we ob- 
tained are 34% and 17% for the flat CMRD and MH09 CMF, 
respectively. The probabilities for all the other samples we 
have considered are given in Table [21 

(jShatskv fc Tokovininl I2002D or with /3=-0.4 for K 
dwarfs primaries (iM azeh et"all 1200311 in the o rbital 
range 0-4 AU. IMetchev fc Hillenbrandl (|2009D and 
IKouwenhoven et al.l (|2005lV already included in the 
sample under discussion, found /3=-0.39±0.36 and 
fi—-Q.2>2> respectively. 

5. DISCUSSION 

The results from the field and Sco 0B2 described in 
section UTT] show an overall trend of CMRDs more peaked 
towards equal mass values than predicted by random 
pairing from Bo2010. This result suggests that the cap- 
ture hypothesis, at least for primary stars in the mass 
range 0.25-7 M© and orbital separation range 1-2400 AU, 
is not the major mechanism for binary formation, as it 
has been alre ady proposed by a large number of previous 
studies (e.g. IClarke fc Pringld ITOQll: IBofRn et al.l UOOSl: 
iBate et al.ll2"003lh 



Fig. 5. — Chi-square best fit Mass ratio distribution for the 
sample of primaries in the field with masses between 0.25-6.5 
M0 over the separation range 1-2400 AU. The best chi-square 
fit is a power-law dN/dq oc g^, with /3 = -0.50 ± 0.29 (x^ = 
0.7 with 7 degrees of freedom). 

Our findings appear to be in agreement with pre- 
dictions from fragmentation theories of binary forma- 
tion, even though we cannot discriminate between dif- 
ferent fragmentation mechanisms. In general, near-equal 
mass binaries are the most likely outcome of fragmen- 
tation in hydrodynamica l simulations (e.g. .Bate ,2000l : 
IKouwenhoven et al.l |2009() . In these simulations, com- 
panions are expected to form by the fragmentation of 
massive accretion regions around stars in ve ry early 
phases of star formation ([Goodwin et al.l [20071 ) and to 
gain mass from the gas reservoir around them. In this 
process, even though a secondary star forms with an ini- 
tial mass close to the opacity limit for fragmentation, it 
will accrete from the circumstellar material, reaching a 
mass roughly similar to the primary (IKouwenhoven et al.l 
2009). These calculations generally predict a relation be- 
tween mass ratio and separation, showing closer binaries 
with higher mass ratios than wider systems (.Bate 2000l 
|2009[ ). This outcome differs from the observational ev- 
idence of a CMRD which is independent of separation 
from few to few thousands AU and suggests that some 
key element is still missing in our models of multiple for- 
mation. 

It should be also noted that the field is likely a mix- 
ture o f systems comi ng from very different environ- 
ments ([Goodwinll2010f) . Binaries might ha ve been pro- 
cessed in different ways ([Parker et al.ll2009l ) and diverse 
star-forming regions may contribute in different degrees 
to the field. Recent res ults from N-body simulations 
([Parker fc Goodwinl[20TI[ l show that the CMRD for very 
low mass binaries is independent of dynamical evolution. 
If this preliminary evidence holds for a broader range 
of primary masses, we can rule out the hypothesis that 
the CMRD was drawn from the IMF at any evolutionary 
stage, suggesting the current CMRD corresponds to the 
birth mass ratio distribution. If this is the case varia- 
tions in the CMRD can be used to trace how different 
star formation regions contribute to the field. 

Interestingly, if we compare the CMRD for solar-type 
primaries in the Pleiades and in the field (MH09) we 



8 



Reggiani et al. 



20 f 



Somple of solor-type stors from Metcfiev & Hlllenbrond (2009) 
Sample of solor-type stars from Krous et al. (2011) 




-0.4 

Log (Mj/M,) 



0.0 



Fig. 6. — Companio n mass ratio distr ibution for solar- 
type stars in Taurus UKraus et al.\ \201i} and in the field 
(MH09). The image shows the comparison between the ob- 
served CMRD in Taurus for the sample of solar type stars 
primaries and the one observed in the field over the common 
range of q. 

obtain a probability of ^1%. Likewise (see also Fig- 
ure [H]), we find a probability of '^1% with the KS test 
between the CMRDs for solar-type primaries in Taurus 
and the field (MH09). Perhaps bound open clusters like 
the Pleiades or extremely low density Taurus-like SFs 
do not con tribute significantly to the field stellar popu- 
lation (e.g. iKroupal Il995t iPortegies ZwartI 120091: lAdamsl 
I2010D . This points toward using binary properties to un- 
derstand which star formation events contribute most to 
the field. Furthermore the comparison between the mass 
ratio distribution in Taurus and Pleiades return a 37% 
KS probability, suggesting they are drawn from the same 
parent population. However the reason why the CMRD 
turns out to be similar in so different environments re- 
mains an open question. 

6. CONCLUSIONS 

We have explored the connection between the CMRD 
and the IMF as a function of primary mass in the field 
and in a few examples of clusters and low density associ- 
ations through Monte Carlo simulations. Using the KS 
test we determined the probability that the two distribu- 
tions are consistent. We have also examined the proba- 
bilities that observed samples are consistent with having 
been drawn from a linearly flat mass ratio distribution 
and a CMF of the form dN/dMz cx M^°-^. Finally we 



have found the best chi-square fit for a composite CMRD 
in the primary mass range 0.25-6.5 M0 {q >0.2). 
Our main results can be summarized as follows: 

- We can reject the hypothesis that the CMRD was 
drawn from the single object IMF for solar type 
stars (28-1590 AU and q > O.I) and M dwarfs (1- 
2400 AU and q >0.2) in the field and A and late-B 
type stars in Sco 0B2 (29-1612 AU and q >0.05). 
The observed CMRDs show a larger number of 
equal-mass systems than would be predicted by the 
IMF. This is in agreement with fragmentation the- 
ories of binary formation. 

- We do not see evidence for variation of the CMRD 
of each sample (M and G stars in the field and AB 
stars in Sco 0B2 association) with orbital separa- 
tion in the ranges explored by the observations. 

- Concerning the observed CMRDs for M dwarfs and 
G stars in the field, we obtain a probability of 6% 
that they are consistent with each other over the 
same range of mass ratios. The CMRD for the 
sample of A and late-B type primaries in Sco 0B2 is 
consistent with both the CMRDs of M and G stars. 
In other words, they all appear to be consistent 
with each other. 

- Regarding the combined CMRD of M and G pri- 
maries in the field and intermediate mass primary 
binaries in Sco 0B2 discussed above over the pri- 
mary mass range 0.25-6.5 Mq, we obtain a chi- 
square best fit following a power law dN/dq oc q^ , 
with (3 — —0.50 ± 0.29, consistent with previous 
studies. 

Certainly further binary studies in young clusters are 
needed to study the dependence of the CMRD on 
dynamical processes and to test possible variations in 
the mass ratio distribution as tracers of different star 
formation mechanisms. 
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